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Abstract. Aeolian sand dunes originate from wind flow and sand bed interactions. According to wind 
properties and sand availability, they can adopt diflterent shapes, ranging from huge motion-less star dunes 
to small and mobile barchan dunes. The latter are crescentic and emerge under a unidirectional wind, with 
a low sand supply. Here, a 3d model for barchan based on existing 2d model is proposed. After describing 
the intrinsic issues of 3d modeling, we show that the deflection of reptating particules due to the shape 
of the dune leads to a lateral sand flux deflection, which takes the mathematical form of a non-linear 
diffusive process. This simple and physically meaningful coupling method is used to understand the shape 
of barchan dunes. 

PACS. 45.70.-n Granular systems - 47.54.-|-r Pattern selection; pattern formation 



1 Characteristics of Barchan Dunes 

R.A. Bagnold opened the way to the physics of dunes with 
his famous book: The physics of blown sand and desert 
dunes From then on, a great deal of investigations - 
laboratory experiments (Ijt2g2j, field measurements j3Sl 
[B1[7||51| ^10II11II12II13II14| and numerical computations |15[ 
Iltilll7lll8lll9ll2(j|l21ll22ll23| - have been conducted by geolo- 
gists and physicists. In particular, a large amount of work 
has been dedicated to the barchan, a dune shaped by the 
erosion of a unidirectional wind on a firm ground. 
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Fig. 1. Side view of a barchan dune . The main properties 
of barchan dunes are outlined: two horns pointing downwind, 
the slip-face and the flat main body. The barchan shown is 
approximately 20 meters long and wide, and 2 meters high. The 
slip-face angle is roughly 30° to the vertical, which corresponds 
to the angle of repose of a sand-pile. 

A side view of a barchan - see Fig. shows a rather 
flat aerodynamic structure. When viewed from above, a 
barchan presents a crescentic shape with two horns point- 
ing downwind. A sharp edge - the brink line - divides the 
dune in two areas: the windward side and the slip face. 



where avalanches develop - see Fig. El Because of a bound- 
ary layer separation along this sharp edge P^E], a large 
eddy develops downwind and wind speed decreases dra- 
matically. Therefore, the incoming blown sand is dropped 
close to the brink line. That is why the barchan is known 
to be a very good sand trapper. Sometimes, when the drift 
of sand is too large, an avalanche occurs and grains are 
moved down the slip face. In short, grains are dragged by 
the wind from the windward side of the dune to the bot- 
tom of the slip-face and, grain after grain, the dune moves. 
Field observations show that barchans can move up to 
70m /year ^1]. Their speed is dependent on wind power 
and on their size: for the same wind strength, the veloc- 
ity of barchans is roughly inversely proportional to their 
heights [T1i^El(Jlll2) . Barchan dimensions range from 1 
to 30 meters high, and from 10 to 300 meters long and 
wide However, large barchans are often unstable, 

leading to complex structures called mega-barchans 
More accurate analyses reveal that the height, width and 
length are related by linear relationships pill3| and that 
no mature barchan dune smaller than one meter high can 
be found: there is a minimal barchan size. Finally, the 
last important characteristic of barchans, but much less 
well-documented, is the sand leak at the tip of the horns 
[Tl | l ||^ 23l. where no recirculation bubble develops. This 
shows that barchan dunes are three dimensional struc- 
tures, whose center part and horns have totally different 
trapping efficiency - see Fig. |21 

Because of the inherent difficulties of field-work, (a typ- 
ical mission duration is rather short compared to the life- 
time of a dune and its shape movement), numerical mod- 
eling is, alongside laboratory experiment Pj, an impor- 
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sand is trapped in the slip lace sand can escape Iroiii Ihe lioms 
Fig. 2. Barchan dune properties. Grains follow the wind direc- 
tion, and sand flux is not strongly deflected by the dune relief. 
As observed in the fleld, sand grains can escape from the horns, 
but not from the main dune body. Instead, they are trapped 
into the slip-face. This difference of behavior between the main 
body and the horns is the key to understand the barchans. 

tant method to explore barchan properties such as their 
morphogenesis, their stabihty and their interaction modes, 
which are still poorly understood. Obviously, these prob- 
lems are related to the original structure of barchans, and, 
accordingly, they should be studied with a 3d approach. 
The aim of the present paper is to discuss how to extend 
existing 2d models to the 3d situation. Then, we will show 
that the crescentic shape can be explained by the exis- 
tence of a sand-flux constraint and a lateral sand-flux. In 
the following part, we start by recalling briefly the main 
features of 2d modeling of dunes . Then a model for the 
3d situation will be presented. 

2 The Cq class of models 

Numerical modeling of dune requires a description of the 
effect of wind flow over a sand-bed. Obviously, comput- 
ing exact turbulent numerical solution starting with the 
Navier-Stokes equations is possible, but this would take a 
very long time. An alternative is to use the model ,23) . 
based on the following approach. 

2.1 Numerical model for the 2d case 

Let us call h{x, t) the sand bed profile and q{x, t) the ver- 
tically integrated volumic sand-flux. They are linked by 
the mass conservation: 

dth + d,q^O. (1) 

The sand flux q{x, t) is the main physical parameter needed 
to understand dune physics. It cannot locally exceed a sat- 
urated value qsat{x,t), which is the maximum number of 
grains that the wind is able to drag per unit of time at x. 
Previous work have already focused on the important role 
played by the saturation of the sand flux and the so-called 
saturation length, I sat Pll22ll23| . The nature of this satu- 
ration process can be understood in terms of sand grain 



inertia: if the wind speed increases, it takes some time for 
the grains, initially at rest, to reach the wind velocity. Dif- 
ferent approaches are possible to describe the evolution of 
the sand flux towards its saturated value qsat- However for 
the sake of simplicity, this saturation process is here taken 
into account with the equation, 

qsat q 

dxq=—j ■ (2) 

This equation keeps the most important effect: the exis- 
tence of a characteristic length-scale Isat- Moreover, qsat, 
depends on the wind shear velocity u* . Even if the nature 
of this dependency is still debated, the saturated flux is a 
growing function of m* (25»,26-27j . Then a linear expansion 
of qsat, as proposed in the innovative work of Sauermann 
et al. |21ll22j from the model of Jackson & Hunt |28j leads 
to: 

^S^ = l + A[^d,K{x-x) + Bd.,K{x), (3) 

Q J X 

where Q is the saturated flux value on a flat ground and 
hf,{x) is the envelope of the dune. This envelope encages 
the dune and its recirculation bubble and is used to 
include the boundary layer separation heuristically. Notice 
that in the latter equation, he{x,t) appears only through 
its flrst spatial derivative. It is consistent with the assump- 
tion that atmospheric turbulent boundary layer is fully 
developed, so changes in u*, and accordingly in qsat are 
scale-invariant. From an aerodynamical point of view, this 
relation takes into account two effects: a pressure effect, 
controlled by A, where the whole shape acts on the wind 
flow; and a destabilizing effect, controled by B, which en- 
sures that the maximum speed of the flow is reached before 
the dune summit. Although in principle we could compute 
the parameters A and B, we prefer to consider them as 
tunable parameters. Finally, avalanches are simply taken 
into account: if the local slope exceeds a critical value, the 
sand flux is increased strongly along the steepest slope. 
We will come back to the description of avalanches in sec- 
tion 4. The only scaling quantities are Q and Isat and they 
are used to adimensionalize the problem. 

In fact, this model belongs to what we called the 
class. It does not depend strongly on the model use for 
the shear stress: other models [2^1 of shear stress pertur- 
bations could be used providing that they include the role 
of the whole shape {A) and the asymmetry of the flow 
[B). The same remarks apply to the charge eauationj21l 
122] and to the equation linking qsat and m*. This shows 
the robustness of the physics ingredients used in this class 
of model. Even though this method does not provide the 
most detailed results, this kind of model has been used 
successfully to model 2d barchan proflle (21..22^23j . 

2.2 Speed dispersion and trapping efficiency 

As a matter of fact, simulations in 2d show the existence 
of two kinds of solutions: dune and dome |2S|- The dune 
solution has a slip-face that catches all the incoming sand: 
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the dune can only grow - except if the input sand flux 
is null. On the contrary, for a dome, a large amount of 
sand can escape, and if the loss of sand is not balanced 
by an influx, the dome can only shrink. Hence, these two 
solutions cannot co-exist. To compute the steady shape of 
a 3d barchan, it is possible to cut it in "slices" parallel to 
the wind direction and, for each slice, to use the 2c? model 
to compute its evolution. However, in a real barchan, a 2d 
slice from the horns seems to behave like 2d dome solution, 
while a 2d slice from the main body coresponds to 2d dune 
solution. However, in order to make slices from the main 
body and slices from the horns coexist, we must introduce 
a sand flux which will redistribuate laterally the sand from 
the center towards the horns. This sand flux coupling is 
also needed to overcome the speed dispersion effect. As 
a matter of fact, if all the slices have initially the same 
shape - but at a different scale ratio - the saturated flux 
at the crest is the same for all the slices, as imposed by 
turbulence scale invariance and the speed of a slice is given 
by: 



where qc and Qout are respectively the flux at the crest 
and the output flux, and he is the height at the crest 
of the slice. Hence, the smaller the slice, the faster its 
motion. This dispersion explains why the barchan takes 
a crescentic shape. But to reach a steady state, all the 
slices must move at the same speed. Therefore, the lateral 
coupling must also induce a speed homogenization. 



3 DifTerent lateral coupling mechanisms 

What physical mechanisms can lead to a redistribution 
of the sand flux on the dune surface? We can think of 
three different possibilities: avalanches, lateral wind shear 
stress perturbations and grain motions (saltation and/or 
reptation) . 

3.1 Avalanches 

First, avalanches develop in three dimensions along the 
steepest slope. This creates lateral sand flux in the slip face 
area. However, for real barchan dunes, that sand-transport 
is directed from the edges towards the center. Accordingly, 
the center part grows and slows down while the border 
slices shrink and accelerates. This does not constitute a 
stabilizing mechanism. 

3.2 Wind deflection 

The lateral wind deflection is another possibility. Field ob- 
servations PP and numerical simulations 12] tend to show 
that, due to its flatness, the dune does not make the wind 
flow deviate too much laterally. Consequently trajectories 
of grains are not dragged into the lateral direction. Nev- 
ertheless, it is always possible to compute the wind speed 



perturbations in the lateral direction and see what it gives. 
This has been done recently '3U' , but there is hardly any 
evidence that it is the main physical process responsible 
for a lateral sand flux on the dune. On the other hand, it 
appears to be an important effect for the development of 
lateral instabilities along transversal dunes [31] ■ 

3.3 Saltation and reptation coupling 



witid direction 




two kind of motions 

incoming salton 




3D case 
lateral deflexion 



Fig. 3. Influence of gravity. While trajectories of saltans axe 
deflected randomly at each collision, the reptons are always 
pushed down along the steepest slope. 

A better candidate is the grain motion. Sand trans- 
portation can be described in terms of two species |27j: 
grains in saltation - the saltans - and grains in reptation 
- the reptons - see Fig. O Saltans are dragged by wind, 
collide with the dune surface, rebound, are accelerated 
again by the airflow, and so forth. At each collision, saltans 
dislodge many reptans, which travel on a short distance, 
rolling down the steepest slope, and then, wait for another 
saltan impact. 

As the wind deflection is weak, we assume that saltans 
follow quasi 2d trajectory in a vertical plane (see last sec- 
tion), except when they collide with the dune. At each 
collision, they can rebound in many directions, depending 
on the local surface properties, and this induces a lateral 
sand flux. Given that the deflection by collision is strongly 
dependent on the surface roughness, we assume that, on 
average, the deflection of saltans is smaller than for rep- 
tans, which are always driven towards the steepest slope. 
In the following derivation, we will neglect the sand flux 
deflection due to saltation collisions. 

Despite the major role played by saltans in dune dy- 
namics, the presence of reptans should not be dismissed. 
According to field observations (32m33) , reptans are strongly 
dependent on the local slope: this can be observed by look- 
ing at the relative orientation of the wind and ripples. On 
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hard ground, ripples are perpendicular to wind direction, 
but on a dune their relative orientations change with the 
local slope |22| ■ Even if one is reluctant to use the saltation 
coupling because of its inherent difficulties, the reptons 
would still be rolling on the dune surface: gravity nat- 
urally tends to move the grains along the steepest slope, 
leading to lateral coupling of the dune slices. In the follow- 
ing part, we will focus on this coupling by reptation, which 
has never been used for the study of barchan dunes. The 
saltation process remains important in this model, since 
it induces reptation coupling. 



4 A model with reptation coupling 

4.1 Formulation of the model 




Fig. 4. Typical initial and final shape of barchan dunes given 
by the 3d Cc model. Parameters used are A — 9.0, B = 5.0, 
D — 0.5 with initially: Wo = SOhat, Ho = Shat- 



Then, the mass conservation equation becomes: 

dth + dxQsal + ^Qrep = 0. (8) 

Calling D = af3/{l + a) and q ~ {l + a)qsah the two latter 
equations can be rewritten as : 

dxQ = —J (9) 

dth + d^q = DidAqdcch) + dy{qdyh)) (10) 

Finally, we obtain the same set of equations as in the 2d 
Cq model, but with one more phenomenological parame- 
ter, D, which can be understood as the importance of the 
lateral coupling, because of lateral deflection of reptons. 
This formulation appears to be a nonlinear difusion equa- 
tion driven by the non dimensionnal coefficient D. For a 
homogeneous flux solution, QD appears to be a diffusion 
coefficient, showing the diffusion-like role of the coupling 
coefficient D. Notice that q is no longer the saltation flux, 
but the part of the flux that does not depend on the bed 
slope. Avalanches are computed in three dimensions us- 
ing a simple trick. If the local slope exceeds the threshold 
fid, the sand flux is strongly increased by adding an extra 
avalanche flux: 

Qa = E{Sfl)Vh, (11) 

where Sfi is null when the slope is lower than fid and equal 
to {6fi ~ |Vft-p — fi^) otherwise. For a sufflciently large 
coefScient E, the slope is relaxed independently of E. Fi- 
nally, note that quasi-periodic boundary conditions (the 
total output ffux is reinjected homogeneously in the nu- 
merical box) are used to perform the numerical simula- 
tions. These boundary conditions are used, flrst to work 
with a constant mass, and second to force the system to 
converge towards its steady state. Obviously these bound- 
ary conditions constrain the final shape. 



As reptons are created by saltons impacts, the part of 
the flux due to reptons is assumed to be proportional to 
the saltation ffux, qsai Hence, the flux of reptons is 
simply written as |34| : 

Qrep = aqsali^x " PV h) (5) 

The a coefficient represents the fraction of the total sand 
flux due to reptation on a fiat bed. In case the bed is not 
flat, the flux is corrected to the flrst order of the h deriva- 
tives, by a coefficient /3 and directed along the steepest 
slope to take into account the defiection of reptons trajec- 
tories by gravity. Assuming that saltation trajectories are 
2d, the total sand fiux q is given by: 

q ^ qsai{i + a)ex - Paqsai'Vh (6) 

Moreover, reptation fiux is assumed to instantaneously fol- 
low the saltation fiux, so that there is no other charge 
equation than the saltation fiux one: 

dqsat qsal /r7\ 
xqsal = -, • (7) 



4.2 The origins of the crescentic shape 

Our 3d Cq model depends on three phenomenological pa- 
rameters : A and B take into account aerodynamics ef- 
fects, and D describes the efficiency of the slices coupling. 
The final shape, if one exists, depends on these parame- 
ters. For example, Fig. 21 shows a typical 3d final steady 
state of a computed barchan dune, which looks likes a real 
aeolian one. The arrows on Fig. [S] indicate the direction 
of the total sand fiux on the whole dune shape. The de- 
viation towards the horns is clearly visible as the sand 
captured by the slip-face. Looking at FigjBlhelps us to un- 
derstand the formation of the crescentic shape of barchan 
dune. Starting from a sand-pile, some horns will appear 
and since they are faster than the center slice of the dune, 
a crescentic shape appears. Now, let us consider a slice 
and let us call qx, the excess flux due to reptons coming 
from the lateral sand fiux, gin the incoming saltation sand 
fiux brought by the wind and qe the fiux due to erosion. 
The existence of a saturated fiux imposes : 

qe + qx + q^n<qTar■ (12) 
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Fig. 6. Barchan formation, a) evolution of an initial cosine bump sand pile : h{x,y) = cos(27ra;/Wo) cos(27rj//W^o)- In the 
beginning, the horns moves faster than the central part of the dune, leading to the formation of the crescentic shape. The shape 
reaches an equilibrium thanks to the lateral sand flux, which feeds the horns. Parameters used are the same than for figures |1] 
and El A = 9.0,B = 5.0,D = 0.5 with initially: Wo = SOhat, Ho = Shat- b) Evolution of a bi -cone sand-pile (same parameters 
than in (a)). This shows that emergence of a crescentic shape is independant on initial shape: the hole in the middle is filled 
up, and a crescentic shape forms. In both boxes the height between two level-lines is 0.2lsat- 
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Fig. 5. Lateral diffusion of sand fiux. The angle of the fiux 
vectors are magnified 3 times to be clearly visible. Both the 
defiection towards the horns and the presence of avalanches 
can be observed. Parameters used are the same than for figure 
Blandlgl A = 9.0,B = 5.0,D = 0.5 with initially: Wo = SOLat, 

Ho = Slsat- 



is governed by the erosion, the slice slows down. After a 
while, the horns receive enough sand from the main body 
to decrease their speed and to compensate the output flux. 
Finally the speed of all slices is the same and the barchan 
moves without changing its shape. Moreover, the center 
part, which would grow without lateral coupling, can now 
have an equilibrium shape, since all the extra flux is devi- 
ated towards the horns. Furthermore this coupling process 
is stabilizing with respect to local deformation. If a slice 
increases in height, the excess sand flux leaving the slice 
will increase, and the deformation will shrink. Similarly, 
starting from a two maxima shape, the part of the flux sen- 
sitive to the local slope tends to fill up the gap between the 
two maxima: the whole mass is redistributed and a single 
barchan shape is finally obtained - see Fig. This agrees 
with the apparent robustness of the crescentic shape ob- 
served on the field. Whatever the external conditions are, 
the same morphology is roughly found everywhere where 
barchans develop. Hence, this lateral coupling helps us to 
understand such barchan structures. However, we have no 
clue about the possible value of this coupling coefficient 
D and it is therefore useful to study the influence of D on 
the barchanic shape. 



where q™"^ is the maximum value of the saturated flux 
on the given slice. Thus, if the excess flux increases, 
the erosion flux, q^, decreases. As the speed of the slice 
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Fig. 7. Different 3d shapes obtained for different values of D 
and for the same initial sand-pile. Parameters used (except for 
D coefficient) are the same than for figures 2] El and |S| Com- 
putations are performed with quasi-periodic boundary condi- 
tions. The barchan central slice is represented on the right. 
Notice that for D = 0, a steady state is reached because of 
lateral sand redistribution due to avalanches on the windward 
side. However, there is a sharp edge which separates the dune 
in two part, leading to an arrow shape rather than to a cres- 
centic shape. 

5 Significance of the coefficients 

5.1 The influence of the coupling coefficient D 

As can be seen on Fig. |3 D has a crucial influence on 
the global morphology. First of all, for a small D, {i.e. a 
small lateral sand flux), the flux escaping from the tip of 
the horns is mainly due to qin and to the erosion flux, 
Qe. Therefore, erosion may be important, and the horns 
may move quickly. Thus, during the transient, the barchan 
elongates. On the contrary, for a larger D, the lateral flux 
received by the horns, qx, is larger and then the erosion is 
less important. The horns move slower and the elongation 
is shorter. 



the width of the horns in the desert may be a convenient 
way to estimate the value of D from field observations. 
Finally the aspect ratio of the dune, h/w, varies with D, 




0.2 0.4 0.6 
Diffusion parameter D 

Fig. 8. Dependance of barchan dune properties on the cou- 
pling coefficient D. The case with D = is taken has a 
reference. The properties are: (a) width of a horn, (b) to- 
tal width, (c) equilibrium input fiux, (d) speed of the dune 
, (e) total length (including horns), (f) length of the center 
slice, (g) maximum height. A typical output fiux cross sec- 
tion is shown in the upper left part of the graph: the width 
of the horns is measured by measuring the width of the fiux 
peaks. The initial sand pile is always a cosine bump defined 
by h{x,y) = Ho coa{2-!Tx/Wo) cos(2-!Ty /Wo) with Wo = 20 and 
Ho = 2.0. 



because the coupling mechanism tends to diminish slopes. 
In other words, the dune tends to spread out laterally 
when D increases, sometimes leading to structures with- 
out a slip-face: 3c? dome solution. In this case, sand can 
escape from the whole transversal section of the dune and 
the output flux is then higher than for the dune solution. 
Basically the transition from 3d dune to 3d dome situa- 
tions occurs when the horn width is equal to half the dune 
width. This corresponds approximately to D ~ 2.0 with 
the parameters A = 9.0 and B = 5.0. For too large a dif- 
fusion, the dune becomes very wide and flat, without a 
slip-face, and remains unsteady. 



Secondly, the width of the horns is also dependent on 
D. For a large D, the deflected flux is large. However, 
as the output sand flux cannot exceed a maximum value, 
horns have to be wide in order to transport away all the in- 
coming sand, and indeed, the width of the horns increases 
with D. This simple explanation is confirmed by numeri- 
cal simulation, as depicted in Fig. Hand Fig.|Hl Measuring 



5.2 The influence of coefficients A and B 

However, looking only at the D influence is not accurate 
enough to describe qualitatively the influence of diffusion. 
In fact, the shape is the result of three physical effects 
governed by the parameters A, B and D. Fig. presents 
a phase diagram with D kept constant at a value of 0.5. 
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Many different morphologies can be observed, from a large 
and thin unsteady crescent to a flat unsteady sand patch, 
and with steady barchans inbetween. First, A governs the 
stability of the sand bed, and then a large A tends to flat- 
ten the sand-pile. Increasing B, on the other hand, forces 
the slope to increase and to nucleate a slip- face. Thus, 
when the ratio A/B is high, there is hardly a slip face. On 
the contrary, when it is low, the slip face develops on the 
whole dune. Thus the A/B parameter controls the way 
the slip-face appears. Secondly, one can observe that the 
dunes obtained with a constant ratio A/B have different 
morphologies. This shows that for small A and B the cou- 
pling coefficient D plays an important role in the shaping 
of the dune. Field measurements of the horns' width and 
the minimal size of barchans should help us to give an es- 
timation of these parameters. The following table provides 
with a summary of the different qualitative effects of A,B 
and D. 



physical mechanism Main effect when increasing 
A curvature effect L increases, 

W decreases. 
B slope effect Slip-face appears easily, 

dxh increases, W increases. 
D Lateral coupling H decreases, 

W and horns width increase, 

output flux increases. 

Finally, the initial sand-pile does not always reach a con- 
stant shape. For large A/B no steady state appears. This 
is due to the fact that the erosion is so strong that it 
cannot be balanced by the influx, and the sand-pile just 
becomes flatter and flatter. On the contrary, for a large 
B/A, a thin unsteady crescent appears. Finally, for a large 
D (compared to A and B) steady and unsteady dome so- 
lutions are produced by the numerical model. This shows 
the importance of A, B and D on the shape and on the 
minimal size of computed barchan dunes. 



6 Conclusion 

We have seen that reptation can be used to describe the 
3d crescentic shape of barchan dunes. As reptons are sen- 
sitive on the local slope, they induce a lateral sand flux. 
This allows a given sand-pile blown by the wind to reach 
an equilibrium shape, despite varying "slices" speed and 
intrinsic differences between the main body and the horns. 
The crescentic shape comes from two mechanisms, which 
compete with each other. The first one is the speed dis- 
persion within the dune: the small height slices are faster 
than higher ones, leading to the crescentic shape, and 
eventually to the destruction of the barchan dune. The 
second one, is the lateral flux deflection that reduces the 
horns' speed by decreasing the erosion flux on them. It 
leads to the homogenization of the speed of the different 
slices, and eventually to the propagation of a steady state 
dune or a steady sand dome. This mechanism is efficient 
since the overall excess flux can escape from the horns. In 
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Fig. 9. Evolution of the same initial cosine bump sand pile 
ih{x,y) = cos(27ra;/Wo)cos(27r2//Wo) with Wo = SOhat, Ho = 
Slaat) changing the value of the parameters A and B. D is kept 
constant and equal to 0.5. The dashed lines separate qualita- 
tively the different stability and shape domains. The dome do- 
main is divided into two parts depending on the steadiness of 
the dome solution. 

other words, barchanic shape appears to be the very basic 
shape taken by any sand-pile blown by a unidirectional 
wind. Now that the steady shape of barchan dunes seems 
to be well described, it would be interesting to study the 
dynamics of such a structure placed far from equilibrium 
conditions, because, in the fleld, input flux has no reason 
to be related to the output flux of a dune. 

As a matter of fact, in this paper, barchans have been 
numerically obtained with quasi-periodic boundary con- 
ditions. If standard boundary conditions are used, a typi- 
cal barchan dune takes the crescentic shape, for the same 
reasons as the ones exposed in this article, but it keeps 
growing or shrinking depending on the incident sand flux, 
either higher or lower than the sand flux escaping from the 
horns 24 . This point is crucial because it suggests that 
single real barchan dunes are not steady structures and, 
moreover, it shows that a crescentic shape can form even if 
the dune is not in an equilibrium state. Therefore, further 
investigations are still required to understand the long- 
term existence and shape of barchan dune in the fleld. 
In our model we have neglected saltation coupling; but 
we believe that it is not the dominant effect. Moreover, it 
could easily be taken into account by adding a new cou- 
pling coefficient dedicated to lateral deflection of saltans. 

Finally, it appears that depending on the coupling co- 
efficient, D, the solution can either be a 3c? dome or a 3d 
barchan. Similarly, changing A and B can lead to different 
shapes, and different behavior, it might be interesting to 
understand the transformation of a 3d sand patch, steady 
or unsteady, into a barchan dune because of fluctuations 
of these parameters, and particularly of D variations. As 
a matter of fact, these parameters may depend on exter- 
nal conditions such as grains size distribution, humidity 
or wind fluctuation. 
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